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Abstract: The gravitino dark matter hypothesis in the brane cosmology is studied. The 
theoretical framework is the CMSSM for particle physics and RS II brane model for gravity. 
It is found that the gravitino can play the role of dark matter in the universe and we 
determine what the gravitino mass should be for different values of the five-dimensional 
Planck mass. An upper bound is obtained for the latter. 
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1. Introduction 

There are good theoretical reasons for which particle physics proposes that new exotic 
particles must exist. Supersymmetry (SUSY) is an ingredient that appears in many theories 
for physics beyond the standard model. SUSY solves the hierarchy problem and predicts 
that every particle we know should be escorted by its superpartner. In order for the 
supersymmetric solution of the hierarchy problem to work, it is necessary that the SUSY 
becomes manifest at relatively low energies, less than a few TeV, and therefore the required 
superpartners must have masses below this scale (for supersymmetry and supergravity see 



One of the theoretical problems in modern cosmology is to understand the nature of 
cold dark matter in the universe. There are good reasons, both observational and theo- 
retical, to suspect that a fraction of 0.22 of the energy density in the universe is in some 
unknown "dark" form. Many lines of reasoning suggest that the dark matter consists of 
some new, as yet undiscovered, massive particle which experiences neither electromagnetic 
nor color interactions. In SUSY models which are realized with R-parity conservation the 
lightest supersymmetric particle (LSP) is stable. A popular cold dark matter candidate 
is the LSP, provided that it is electrically and color neutral. Certainly the most theoret- 
ically developed LSP is the lightest neutralino |||. However, there are other dark matter 
candidates as well, for example the gravitino |3|, |[] and the axino ||, the superpartner of 
axion || which solves the QCD problem via the Peccei-Quinn mechanism |?J|. In this article 
we work in the framework of Randall-Sundrum type II brane model (RSII), we assume that 
the gravitino is the LSP and address the question whether the gravitino can play the role 
of dark matter in the universe, and for which range for gravitino mass and five-dimensional 
Planck mass. We remark in passing that gravitino production in brane cosmology has also 
been considered in Q. 



e.g. 1). 
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For a better qualitative understanding of the present work, we report here the main 
outcome of this article. The usual gravitino overproduction problem can be avoided because 
of the faster cosmological expansion in early RSII cosmology. Thus, the gravitino becomes 
a viable dark matter candidate for a wide range of parameters of the model. This requires, 
in particular, that the fundamental Planck mass M5 is not too high. A key quantity in our 
analysis is the transition temperature Tt, which is determined solely by M5. Both thermal 
and non-thermal gravitino production mechanisms are affected by the novel expansion law 
of the universe and are modified as follows. In standard four-dimensional cosmology the 
gravitino yield from thermal production is proportional to the reheating temperature Tr 
after inflation. However, in RSII brane model considered here the gravitino yield becomes 
proportional to the transition temperature. Furthermore, the gravitino yield from non- 
thermal production differs from the standard expression by a factor which is inversely 
proportional to T t . 

Our work is organized as follows: The article consists of four section, of which this 
introduction is the first. In the second section we present the theoretical framework, while 
in section 3 we show the results of our analysis. Finally we conclude in the last section. 



2. The theoretical framework 



2.1 The brane model 



Over the last years the brane- world models have been attracting a lot of attention as a novel 
higher-dimensional theory. Brane models are inspired from M/string theory and although 
they are not yet derivable from the fundamental theory, at least they contain the basic 
ingredients, like extra dimensions, higher-dimensional objects (branes), higher-curvature 
corrections to gravity (Gauss-Bonnet) etc. Since string theory claims to give us a funda- 
mental description of nature it is important to study what kind of cosmology it predicts. 
Furthermore, despite the fact that supersymmetric dark matter has been analyzed in stan- 
dard four-dimensional cosmology, it is challenging to discuss it in alternative gravitational 
theories as well. Neutralino dark matter in brane cosmology has been studied in ||, while 
axino dark matter in brane- world cosmology has been studied in [10]. 

In brane-world models it is assumed that the standard model particles are confined 
on a 3-brane while gravity resides in the whole higher dimensional spacetime. The model 



first proposed by Randall and Sundrum (RSII) [11|, is a simple and interesting one, and 
its cosmological evolutions have been intensively investigated. An incomplete list can be 
seen e.g. in [O]. In the present work we would like to study gravitino dark matter in the 
framework of RSII model. According to that model, our 4-dimensional universe is realized 
on the 3-brane with a positive tension located at the UV boundary of 5-dimensional AdS 
spacetime. In the bulk there is just a cosmological constant A5, whereas on the brane there 
is matter with energy- momentum tensor t^ iu . Also, the five dimensional Planck mass is 
denoted by M5 and the brane tension is denoted by T. 

If Einstein's equations hold in the five dimensional bulk, then it has been shown in 
1 13] that the effective four-dimensional Einstein's equations induced on the brane can be 
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written as 



+ A 4 5 Ml/ = — 2"T^ + (773) T^iw - E^y (2.1) 

where 5^ is the induced metric on the brane, -k^ = T2 t t ^v + \ 9^v T a /3 T a ^ — \ T/mx t v ~ 
21 1-2 9p.v 1 A4 is the effective four-dimensional cosmological constant, m p i is the usual four- 
dimensional Planck mass and = C^ pa n a n p g^ g" is a projection of the five- dimensional 
Weyl tensor C a p pa , where n a is the unit vector normal to the brane. The tensors TVu U and 
describe the influence of the bulk in brane dynamics. The five-dimensional quantities 
are related to the corresponding four-dimensional ones through the relations 




m pl = A\l— Ml (2.2) 
and 

A *=2lf( A5 + ^f) (2 ' 3) 

In a cosmological model in which the induced metric on the brane g pu has the form of a 
spatially flat Friedmann- Robertson- Walker model, with scale factor a(t), the Friedmann- 
like equation on the brane has the generalized form [jO ] 

TT o A4 8^ 1 n C . 

where C is an integration constant arising from £w. The cosmological constant term and 
the term linear in p are familiar from the four-dimensional conventional cosmology. The 
extra terms, i.e the "dark radiation" term and the term quadratic in p, are there because 
of the presence of the extra dimension. Adopting the Randall-Sundrum fine-tuning 

T 2 

A 5 = o (2.5) 

5 6M| 

the four-dimensional cosmological constant vanishes. In addition, the dark radiation term 
is severely constrained by the success of the Big Bang Nucleosynthesis (BBN), since the 
term behaves like an additional radiation at the BBN era |p^| . So, for simplicity, we neglect 
the term in the following analysis. The five-dimensional Planck mass is also constrained 



by the BBN, which is roughly estimated as M5 > WTeV [15]. The generalized Friedmann 
equation takes the final form 



r2 %kGJ, , P 



5=^1 + - (2-6) 
3 V PoJ 

where 

po = 96nGMl (2.7) 

with G the Newton's constant. One can see that the evolution of the early universe can be 
divided into two eras. In the low-energy regime p <C po the first term dominates and we 
recover the usual Friedmann equation of the conventional four-dimensional cosmology. In 
the high-energy regime po <C p the second term dominates and we get an unconventional 
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expansion law for the universe. In between there is a transition temperature Tt for which 
p(Tt) = pq. Once M5 is given, the transition temperature Tt is determined as 



T t = 1.6 x W< t1 GeU (2.8) 

where <7 e // counts the total number of relativistic degrees of freedom. 
2.2 The particle physics model 

The extension of standard model (SM) of particle physics based on SUSY is the minimal 
supersymmetric standard model (MSSM) fl6[| . It is a supersymmetric gauge theory based 
on the SM gauge group with the usual representations (singlets, doublets, triplets) and 
on M = 1 SUSY. Excluding gravity, the massless representations of the SUSY algebra are 
a chiral and a vector supermultiplet. The gauge bosons and the gauginos are members 
of the vector supermultiplet, while the matter fields (quarks, leptons, Higgs) and their 
superpartners are members of the chiral supermultiplet. The Higgs sector in the MSSM 
is enhanced compared to the SM case. There are now two Higgs doublets, H u ,Hd-, for 
anomaly cancelation requirement and for giving masses to both up and down quarks. After 
electroweak symmetry breaking we are left with five physical Higgs bosons, two charged 
and three neutral A, H, h (h being the lightest). Since we have not seen any superpartners 
yet SUSY has to be broken. In MSSM, SUSY is softly broken by adding to the Lagrangian 
terms of the form 

• Mass terms for the gauginos gi, M\, M3 

Mgg (2.9) 



Mass terms for sfermions / 



m)pf (2.10) 



• Masses and bilinear terms for the Higgs bosons H u ,H d 

m 2 Hu HlH u + m 2 Hd H ] d H d + Bp{H u H d + h.c.) (2.11) 

• Trilinear couplings between sfermions and Higgs bosons 

AYf\Hf 2 (2.12) 

In the unconstrained MSSM there is a huge number of unknown parameters [17] and thus 



little predictive power. However, the Constrained MSSM (CMSSM) or mSUGRA (!§] is 
a framework with a small controllable number of parameters, and thus with much more 
predictive power. In the CMSSM there are four parameters, mo, mi/21 A)> tanf3, which are 
explained below, plus the sign of the fi parameter from the Higgs sector. The magnitude of 
[i is determined by the requirement for a proper electroweak symmetry breaking, its sign 
however remains undetermined. We now give the explanation for the other four parameters 
of the CMSSM 
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• Universal gaugino masses 



M^Mgut) = M 2 (M GUT ) = M 3 (M GUT ) = mi/2 



(2.13) 



Universal scalar masses 

m^(M GUT ) = m 

Universal trilinear couplings 

A^Mgut) = Afj (Mqut) = Aj(Mgut) = 4><*v 



(2.14) 



(2.15) 



tan/3 



v 2 



where v\ , v 2 are the vevs of the Higgs doublets and Mqut 
Unification scale. 



(2.16) 

10 16 GeV is the Grand 



3. Analysis and results 

We consider four benchmark models (shown in the table below) for natural values of 
7710,7711/2, representative values of tan(3 and fixed A$ = 0, fi > 0. In these models the 
lightest neutralino (denoted by x) is the lightest of the usual superpartners and thus the 
NLSP. Furthermore the following experimental constraints (for Higgs, scalar r and chargino 
masses and a rare decay) [0, 19 1 are satisfied 



m h > 114.4 GeV 



m fl > 81.9 GeV 
m-± > 81.9 GeV 

Al 



BR(b -» s 7 ) 



(3.39i°;f 7 ) x 10- 4 



(3-1) 
(3-2) 
(3-3) 
(3-4) 



Model 


m (GeV) 


"2-1/2 (GeV) 


tan(3 


m x (GeV) 


n x h 2 


A 


200 


500 


15 


205.42 


0.64 


B 


400 


800 


25 


337.95 


1.82 


C 


1000 


600 


30 


252.41 


7.37 


D 


350 


450 


20 


184.46 


1.2 



Table 1: The four benchmark models considered in the analysis 



The SUSY spectrum (as well as the Higgs bosons masses) is computed using the 
FORTRAN code SuSpect |2(| , the neutralino relic density is computed using the computer 
program micrOMEGAs [^l[ and the top quark mass is fixed to mt = 172.7 GeV | p2]| . 
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For the gravitino abundance we take into account both thermal (TP) and non-thermal 
production (NTP) 

n 3/2 h 2 = n^ p h 2 + n^f 2 h 2 (3.5) 

and we impose the WMAP constraint for cold dark matter [^] 

0.075 < n cdm h 2 = n 3/2 h 2 < 0.126 (3.6) 

In the NTP case the contribution to the gravitino abundance comes from the decay of the 
NLSP 

^3/2^ 2 = ^NLSph 2 (3.7) 

m NLS p 

with m.3/2 the gravitino mass, tunlsp the mass of the NLSP and &NLSph 2 the NLSP 
abundance had it did not decay into the gravitino. In the TP case the gravitino abundance 
is given by 

where To = 2.7 K is the temperature of the universe today, s = /i e //27r 2 /45 is the entropy 
density for the relativistic degrees of freedom, p cr = (3Hq)/(8ttG) is the critical energy 
density today and Y 3 / 2 is the gravitino yield defined by Y 3 / 2 = n 3 / 2 /s, where n 3 j 2 is the 
gravitino number density. For the yield in standard four-dimensional cosmology we employ 
the formula (3) of Q, which updates |2jj ] an older calculation p5fl . We see that y^f(To) 
is proportional to the reheating temperature after inflation. 

Now we need to take into account the effects of the novel law for expansion of the 
universe. These are twofold: First, in the TP case the reheating temperature is replaced 
by the transition temperature, Tr — ► 2T t [10, 26]. Second, the relic density of a particle of 



mass m is modified as follows [27] 



n&- c J^ (3 ' 9) 

where C is a numerical factor of order one, the index b stands for "brane", the index s 
stands for "standard", xt = m/Tt and = m/Td, with Tt the transition temperature and 

(s) 

Td the decoupling temperature of the particle of mass m. In standard cosmology x d ~ 30. 

For each benchmark model we have obtained plots (below we show as an example 
the plots for model A, for the rest of the models there are similar plots) which show the 
gravitino abundance fl 3 / 2 h 2 as a function of the gravitino mass m 3 / 2 for several values of 
the five-dimensional Planck mass M5. We discriminate between low values and high values 
of M5. In the TP case, the calculation assumes that the particles are in thermal equilibrium 
and thus the analysis is valid for temperatures T > 1 TeV or for five-dimensional Planck 
mass M5 > 1.83 x 10 8 GeV. We use this criterion to divide the values of M5 into two 
categories, low values M5 < 1.83 x 10 8 GeV and high values M5 > 1.83 x 10 8 GeV. Figure 
1 corresponds to the first case (the values that we have used are 5 x 10 GeV, 8 x 10 5 GeV 
and 6 x 10 6 GeV from top to bottom), while figure 2 corresponds to the second case (the 
values that we have used are 10 10 GeV, 10 11 GeV and 10 12 GeV from bottom to top). In 
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the low M5 case there is always one allowed range for the gravitino mass from the sub-Gey 
range to a few tens of GeV. On the other hand, in the high M5 case there are two allowed 
ranges for the gravitino mass, one in the sub-Gey range and another one in the -a few 
GeV up to a few tens of GeV- range. If however M5 is high enough, M5 > 10 12 GeV, there 
is no allowed range at all. 

4. Conclusions 

We have studied gravitino dark matter in the brane-world cosmology. The theoretical 
framework for our work is the CMSSM for particle physics and RS II for gravity, which 
predicts a generalized Friedmann-like equation for the evolution of the universe. We as- 
sume that gravitino is the LSP and the lightest neutralino is the NLSP. For the gravitino 
abundance we have taken into account both thermal and non-thermal production and have 
imposed the cold dark matter constraint 0.075 < VL c d m h 2 < 0.126. The formulae valid in 
standard four-dimensional cosmology are corrected taking into account the novel expan- 
sion law for the universe. We have considered four benchmark models for natural values 
of mo and m 1( / 2 and representative values of tan(3. In these models the neutralino is the 
lightest of the usual superpartners (and thus the NLSP, since we assume that the gravitino 
is the LSP) and experimental constraints are satisfied. For each benchmark model we have 
produced plots of the gravitino abundance as a function of the gravitino mass for several 
different values of the five-dimensional Planck mass. The obtained plots show that in gen- 
eral the gravitino can be the cold dark matter in the universe for gravitino masses from a 
few GeV up to tens of GeV. Furthermore an upper bound on the five-dimensional Planck 
mass is obtained, M 5 < 10 12 GeV. 
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Figure 1: Gravitino abundance versus gravitino mass (in GeV) for several values of the five- 
dimensional Planck mass (benchmark model A and low M 5 ). The strip around 0.1 is the allowed 
range for cold dark matter. Values of M 5 used arc 5 x 10 4 GeV, 8 x 10 5 GeV and 6 x 10 6 GeV 
from top to bottom. 
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Figure 2: Same as figure 1, but for high M5. Values of M5 used are 10 10 GeV, 10 11 GeV and 
10 12 GeV from bottom to top. 
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